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The dynamics of fission has been formulated by generalising the asymptotic expansion of the 
Fokker-PIanck equation in terms of the strength of the fluctuations where the diffusion coefficients 
depend on the stochastic variables explicitly. The prescission neutron multiplicities and mean ki- 
netic energies of the evaporated neutrons have been calculated and compared with the respective 
experimental data over a wide range of excitation energy and compound nuclear mass. The mean 
and the variance of the total kinetic energies of the fission fragments have been calculated and 
compared with the experimental values. 
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^ , I. INTRODUCTION 

At present, it is commonly agreed upon that the fission process is a dissipative phenomena, where initial energy of 
■ the collective variables get dissipated into the internal degrees of freedom of nuclear fluid giving rise to the increase in 
, internal excitation energy. As dissipation is referred to the interaction of the system coordinate with the large number 
of degrees of freedom of the surrounding reservoir, this process is always associated with the fluctuations of relevant 
" ' ' physical observables. Thus, the dynamics of fission process resembles the standard Brownian motion problem, where 
^ I ' the collective variables such as shape degrees of freedom act as 'Brownian particles' interacting stochastically with 
[ large number of internal nuclconic degrees of freedom constituting the surrounding 'bath'. This mesoscopic description 
is inevitable once the fluctuations of the observables are amenable to experimental observation. 
T-H There have been several attempts in the past to study the dynamics of fission by solving either the Langevin 

, equation or multidimensional Fokker-PIanck equation ||^^, which is a differential version of Langevin equation. 

In the case of fission, it is experimentally observed that the variances of the physical observables are, in general, 
small compared to their respective mean values ( typically, the ratio of the root mean square deviation and the 
mean of the kinetic energy is 0.1). The question naturally arises whether one can utilise this simple fact in the 
theoretical scheme instead of solving the Langevin equation (LE) or corresponding Fokker-PIanck equation in detail. 
In this spirit, we present an alternative theoretical prescription for the calculation of various moments of the physical 
^ ' observables related to the fission process based on the assumption that the full solution of the Fokker-PIanck equation 
(~| (FPE) admits an asymptotic expansion in terms of strength of the fluctuations. The asymptotic expansion method 
^ was first developed by van Kampen ^ for the stochastic processes having constant diffusion coefficients. However, 
• w^ , a generalisation of the above prescription is necessary in the case of fission where the dissipation is usually assumed 
' to depend on the instantaneous shape of the fissioning system and therefore the diffusion coefficients are also shape 
5h \ dependent. To the best of our knowledge, such an application in the case of fission is not available in the literature. 
. . 1 In the present paper, we report a generalised formulation of the asymptotic expansion of the Fokker-PIanck equation 
where the diffusion coefficient depends on the stochastic variables explicitly. In this formulation the dynamics of the 
stochastic processes reduces to a set of linear ordinary differential equations which are far simpler to solve as com- 
pared to either multidimensional partial differential ( Fokker-PIanck) equations or stochastic differential (Langevin) 
equations. Presently, we apply this formulation to calculate the various moments of the relevent physical observables 
of the fission process. 

The paper is organised as follows. The generalised formulation and its application to the dynamics of fission is given 
in Sec. II. The calculations and numerical results are discussed in Sec. III. Finally, concluding remarks are given in 
Sec.IV. 

II. ASYMPTOTIC EXPANSION OF THE FOKKER-PLANCK EQUATION 

A. The Formalism 



o 



The mesoscopic description of the fission process begins with a set of Langevin equations: 
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X, = h,{{X},{Y})+r],{t) (2.1a) 
Y,=H,{{X},{Y}) ;z^l,...,N (2.1b) 

where hi and Hi are given functions of the stochastic collective variables Xi, X2, Xn and Yi, I2, in the 

fission process and rji (t) refers to the driving noise term associated with the interaction of the ith collective variable 
with the reservoir constituting nucleonic degrees of freedom. For simplicity, we assume the noise to be a gaussian 
white with zero mean and decoupled for different degrees of freedom with auto-correlation functions given by 

< n,{t) >= 0, < V^{t)V3{t') >= D,{yC}5{t - t')6,j, (2.2) 

where Di{yi) is the diffusion coefficient associated with ith variable, depending only on the sample space yi for the 
stochastic variable Yi. 

The Fokker- Planck equation corresponding to the Langevin equation (2.1) is 

d.n{x},{y},t) ^, d{hj) , djHJ) d^f 

di - + - (l/2)A(y.)^]. (2.3) 

The quantity fi{x},{y},t) is the probability density function depending on the variables xi, X2, s^Af, J/i, 2/2, 
and time t explicitly. If we are interested in finding the time evolution of the conditional probability distribution 
function then we have to solve Eq.(2.3) with initial values Xi{0) = x^, yi{0) = y^, Vi, at t = 0. That is, we have to 
solve Eq.(2.3) for those realisations which are known to start from these specific points in the whole sample space. 

In the cases where diffusion coefficient is constant the asymptotic expansion method of van Kampen consists of 
writing the stochastic variables as the sum of deterministic value and a fluctuating part at each time t with root of 
the diffusion constant as a strength of the fluctuating part. In the present paper, we generalise this method for the 
situations where the diffusion coefficients depend on the stochastic variables explicitly. Such a situation is encountered 
in the case of fission process, where the friction coefficient depends explicitly on the collective variable or shape of the 
nucleus at each instant of time. In this case, we further assume that, in the asymptotic expansion, the strengths of 
the fluctuating parts of the stochastic variables depend only on the deterministic values of the respective y variables: 

Xi = Xi + y/D{y,)Q (2.4a) 

y^ = y^ + VD{y^)^^ (2.4b) 

The quantities {Ci}, {^i} refer to the fluctuations of the stochastic variables {xi} and {yi} around their deterministic 
values {xi},{yi}. Next, we introduce the new distribution function Q depending only on the variables {Cj}i{^j} ^^^id 
t. The normalisation condition suggests that the / and Q will be related by 

N 

/({x}, {y}, t) = l[[DM)r'Q,{C^,^^,t) (2.5) 

i=l 

Substituting Eq.(2.4) in the Fokker-Planck equation(2.3), making Taylor expansion of h{{x},{y}),H{{x},{y}) 
around {a;}, {y} and collecting coefficients of various order of D{yi),we could generate a hierarchy of equations. As 
expected, the first set would give rise to the equation of motion for {x} and {y}. 

Xi = hi{{x} , {y}) (2.6a) 
y\=H,{{x},{y}) -yi (2.6b) 

Eqs.(2.6) are the Euler-Lagrange equation for deterministic motion. These equations are to be solved with initial 
conditions {a;(0)} = {x°},{y(0)} = {y"}. Next, we are going to calculate the conditional probability distribution 
fi{x}, {y}, t I {x"}, {yO}, 0) or Q({C}, {O, ^ I 0, 0, 0). 

Assuming the variation of diffusion coefficient over the narrow width of the distribution function at any instant 
of time to be 0(^)i we could replace the second Fokker-Planck coefficient D{y) by D{y) at each instant of time. 
This assumption makes the calculation extremely simple. Collecting coefficients 0(^°)j we get back quasilinear 
Fokker-Planck equation for Q: 

where ai,bi,Ci, di are given by 
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dh Djyi) 

6. = (^) (2.8b) 
. = S-(M) (2.80 



, dH 



di = (^) (2.8d) 



Eq.(2.7) suggests that 



Qm,{0,t) = l[Qj{Cj,^j,t) (2.9) 
j 

where the distribution function Qj for each j satisfies the similar equation written below without the subscript: 

^ + (M)n - _r«^ + 5^ + + - (1/2)^1 (2 10) 
subject to the initial condition 

Q(c,e,i = o) = 5(0^(0 (2.11) 

The solution of Eq.(2.10) is given by 

QiU.t)'[j^\j /e.pM.K + .« + fe^2^!±£W + 2CWM,,M, (2.2, 
where g{t), G{t), C{t) satisfy the set of coupled first order differential equations : 



with the initial conditions 



5(0) = G(0) = C(0) = (2.14) 



Once Q{C,,£,:t) is known, from Eq.(2.9) and Eq.(2.5) the full conditional probability distribution function 
/({x}, {y}, t I {y"}, 0) is known. Integrating this function over all variables except one, say Xi, one identi- 

fies gi{t) as the variance of the stochastic variable Xj. 

<{Xi-<Xi>f>=gi{t) (2.15) 

Wc note that the homogeniety of Eq.(2.10) suggests that < (^{t) >=< ^{t) >= 0, or the average of the variables 
X and Y at any time will be determined by their deterministic values obtained by solving Euler-Lagrange equation 
(2.6). Similarly, one observes from Eq.(2.12), 

< {Xi- < Xi >){Yi- <Yi>)> = Ci{t) (2.16a) 
< {Yi- < Yi >f > = Gi{t) (2.16b) 
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B. Application to the fission process 
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In the fission process, in accordance with our previous work pO| , the shape of the fissioning nucleus is described 
in terms of the elongation axis (the neck parameter of taken equal to zero). Thus, in the dynamical description 
we have the elongation axis, its relative orientation with respect to an inertial system and respective velocities 
associated with them as the stochastic variables interacting with a large number of internal nucleonic degrees of 
freedom constituting a heat bath at temperature T determined by the excitation energy available to it. We further 
assume that the 'collisional' time scale of the nucleonic degrees of freedom is much shorter than the time scale of 
the macroscopic evolution of the collective variable so that at each instant of time the heat bath is assumed to be in 
quasi-stationary equilibrium. 

The Euler-Lagrange equations (2.6) were solved in our earlier works [|l0| . To avoid repetition we deliberately omit 
the procedure and scheme to solve those equations. For the sake of completeness we merely write those equations and 
refer to our previous papers to clarify the details. 

Giving correspondence to the terminology used in this paper, we associate 

Y^=r,Xi=f, (2.17a) 
Y2^e,X2^e. (2.17b) 

Thus we have 

(2.18a) 

(2.18b) 

(2.18c) 
(2.18d) 
(2.18e) 
(2.18f) 

The quantities Vc, Vn represent the Coulomb and nuclear interaction potentials and 7, 74 are the radial and 
tangential components of friction, respectively. The nuclear part of the interaction is approximated by the proximity 
interaction [ p^ . /i, I2 are the moments of inertia of the two lobes and L refers to the relative angular momentum, gi 
and (72 are the distances of the centres of mass of the two lobes from the centre of mass of the composite dinuclear 
system and the term [172(^2 — 0)+ 9i{0i — 6)] represents the relative tangential velocity of the two lobes. The quantities 
/ and ^ are the moment of inertia and the reduced mass associated with the fissioning liquid drop, respectively ]l3t . 

It has already been shown that the tangential friction which causes dissipation of relative angular momentum L 
into the angular momenta Ii , I2 of the two fragments does not have any significant effect on the physical observables 
. Besides, no experimental observation of angular momentum dispersion of fission fragments are available in 
the litterature. Therefore, in the following the calculations of higher moments are restricted to the radial degree of 
freedom only. For the sake of convenience we omit the subscripts in the functions H and h below. 

The variances are obtained by solving Eqs.(2.13). The diffusion coefficient D is evaluated employing Einstein's 
fluctuation dissipation theorem. Thus, the Eqs.(2.13) now becomes, 

m = 2(|^)g(t) + 2(|^)C(i) + 2j{r)T{r)/fi' (2.19a) 

G{t) = 2C{t) (2.19b) 
rlh r)h 

C{i) = {^)C{t) + 2{^)G{t) + g{t) (2.19c) 
with the initial conditions(2.14). The initial conditions of r and r for solving Eq.(2.18) are ||l^ 

ro = r{t = 0) = r„,„ ± <5ro, ro = r(t = 0) = [^^f/^ (2.20) 

where the potential energy surface around the minimum is approximated as harmonic oscillator with oj being the 
oscillator frequency. At each instant of time, it is assumed that the state of the nucleus is ameanable to a thermody- 
namic description with temperature T. Therefore, 5ro in Eqn. (2.20), which is taken as the root of the thermal average 
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of mean quantum dispersion around the minimum of the potential, is expressed as Sro = "^o^^^ 2T(r^'r ■ ) )"^^^- 
The quantity i?Ar is a random number between and 1 from uniform probabihty distribution and Eq is the initial 
available energy. The temperature T(r) in Eqn. (2.19a) has been calculated from the instantaneous excitation energy 
E*{r) using the relation T(r) = y/{E*{r)/a) with a = A/10. At each instant, the dissipated energy is added to, and 
the energy carried away by the prescission particles (if any) is subtracted from, the excitation energy which is then 
used to calculate the temperature at the next instant. 

Solving Eq.(2.18) and Eqs.(2.19) simultaneously with initial conditions (2.14) and (2.20) we generate the conditional 
probability distribution function /(r, f,i | r{t — 0),r(t — 0),0). The probability distribution function f(r^f,t) could 
be obtained as 

/(r, f,t)= J J{r, f, t\r(t^ 0), r(t = 0), Q)f(r{t = 0), f(t = 0), Q)dr{t = {))df{t = 0) (2.21) 

where f{r{t — 0),r(t = 0),0) is the probability distribution of position and velocity of the stochastic variables at 
the initial time. As described by the initial condition(2.20), this can be represented as 

f{r{t = 0), f(t = 0), 0) = 6{r{t = 0) - r™„ T ^o) x f{f{t = 0)) (2.22) 

Here, we assumed that each fissioning nucleus in the ensemble starts from a fixed initial position but with different 
partioning of initial excitation energy ||l^. Finally, substitution of Eq.(2.22) in Eq.(2.21) would give 

/(r, r, i) = E ^ I (^^^)'^'' 0) (2.23) 

The asymptotic expansion in our model thus provides the following picture: In the zeroth order approximation, the 
motion is described by the Euler - Lagrange equation. This requires the initial momentum as a generator of motion, 
which is supplied by a random fraction of initial available energy Eq of the system. This initial randomness restricts 
the trajectories to have fission fate thus providing the crosssection of the residue. The first order approximation 
provides mostly the other transport property, namely the variance of the physical variable. As the approximation 
of the distribution function is over the solution of the Euler - Lagrange equation of motion, this part provides the 
observables of the escape part of the distribution. In this way, this model could describe the bifurcation of the total 
distribution function if one would solve the full Fokker - Planck equation or the Langevin equation. 



III. NUMERICAL CALCULATION AND RESULTS 

The applicability of the generalised formalism developed in Sec. || has been tested quite rigorously by confronting 
it with a wide range of experimental data on various physical observables of the fission process. The details of such 
calculation procedure has been reported elsewhere (Ref. JlO|), and is given here in brief. 



A. The shape, friction and dynamics of the fissioning system 

For the present calculation, instantaneous shape of the fissioning nucleus is taken to be of the form 10 lUl, 



p\z) = c-\c^ - 7J){A + B7J + azc), (3.1) 

where the coefficients A and B are given by, A = — and B = {c — l)/2, respectively. The variable c 

corresponds to the elongation and a is a parameter which depends upon the asymmetry {aasy) defined as aasy = 
[Ai — A2)/AcN, where Acn is the compound nucleus mass, and Ai , A2 correspond to the masses of the two fragments. 
The parameter a is related to the asymmetry aasy through the relation a = .11937a^sj^ + .24720aQsy jlOtl. As the shape 
changes gradually, the coordinates of the two maxima and that of the minimum of the surface Eq. ( |3.1| ) change. The 
scission point is defined when the minimum point touches the z-axis and it is given hy A — c^a^ /AB — . Therefore, 
the value of c at which scission occurs depends on a and the dependence is given by Csc = — 2.0a^ + .032a + 2.0917. 

The variable r is defined as the centre to centre distance between the two lobes. From the generalised shape given 



by Eqn. (3.1), we first construct the centres of mass of left and right lobes, and call them zi and Zr respectively. Then 



r is defined as r = \zi — Zr\ ■ The reduced mass parameter /x, is obtained from the calculated masses of the two lobes. 
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The temporal evolution of shape of the fissioning nucleus is assumed to start from the minimum of the potential 
energy surface eventually leading to scission. The fission trajectories are obtained by solving the Euler-Lagrange equa- 
tions with conservative forces derived from the nuclear and Coulomb potentials For the non-conservative 
part of the interaction, we would consider viscous drag arising not only due to two body collision but also due to the 
collisions of the nucleons with the wall or surface of the nucleus. Hence 7 in Eqn.(2.18b) contains two parts; 7^^ and 
7*^^, for two-body and one-body dissipative mechanisms, respectively. Assuming the nucleus as an incompressible 
viscous fluid, and for nearly irrotational hydrodynamical flow, 7^^ is calculated by use of the Werner- Wheeler method 
and is given by 

7^^ = RcNh{a)f{ — ) / dzp2[3^^2 ^ p^A^^] (3.2a) 



where the factor h(a) = exp(— iiTa^) is included in order to explain the observed fragment asymmetry dependence 
of neutron multiplicity (for details, see Ref. fl^), and, 

M^) = -^^§-Jj^W). (3.2b) 

The quantities A'^,A'^ are the first and second derivatives of Ac{z) with respect to z. /io is the two body viscosity 
coefficient. The factor /(ff ) is taken to be 

where x — r/RcN, Rcn being the radius of the compound nucleus. 
The tangential friction 7^^ is calculated using the following relation ||l0[| , 

ir = i^r^r, (3.4) 

where n (the value of p at the minima of in Eqn. ^^) is the instantaneous neck radius of the fissioning system. 
One body dissipative force, Fdis, is obtained from the rate of energy dissipation, Edis, by 

Fd^s = -^Ems{x) (3.5) 

where x refers to the rate of change of x with respect to time and Edis{x) is the rate of energy dissipation at x 
given by 

1 _ / .2 

Edis = -^PmV (b dS en , (3.6) 

where e„ is the unit normal direction at the surface. The integration is done over the whole surface, pm is the 
nuclear density and v is the average nucleonic speed obtained from the formula 

v = ^{ — ){EaJaf'^ (3.7) 

with Eg^ is the available energy and the level density parameter, a, is taken to be ^cjv/lO. For the generalised 



shape (3T), one-body friction, 7*^^, is obtained as 

7°^ = 27r pmvRlj,f{dc/dx) / dz p [I + pi\^l\A,pi + (1/2) pPl^f (3.8) 



where p/, are the derivatives of p, A^ with respect to z and all other quantities are defined earlier. The tangential 



part of the one-body friction is calculated in a similar manner as in Eqn. 3.4. 

The friction forces used in the calculation are taken as follows One-body 'wall' friction has been used in the 
ground state to saddle region, where nuclear shapes are nearly mononuclear. The strength of the one-body friction 
used was attenuated to 10% of the original 'wall' value. This weakening of the wall friction has also been confirmed 
from the study of the role of chaos in dissipative nuclear dynamics jlTj l. In the saddle to scission region, on the 
other hand, the nuclear dissipation was taken to be of two-body origin and the value of the viscosity coefficient 
po used in the present calculation was (4 xlO~^^MeT^ . sec • fm~^). This value of pa corresponds to 0.06 TP 
(ITP = 6.24 X IQ-^^MeV ■ sec ■ Jm,-^). 
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B. Prescission neutron emission 



1. Prescission Neutron Multiplicities 

The emission of the prescission neutrons is simulated in the following way. During the temporal evolution of the 
fission trajectory the intrinsic excitation of the system, and vis-a-vis, the neutron decay width at each instant, r„, is 
calculated using the relation r„ = hWn- The decay rate Wn is given by, 

where, (PUn/dEdt is the rate of decay A ^ A— l + ninan energy interval [E,E-{-dE] and a time interval [t,t + dt]. 
The quantity d^Iln/dEdt may be evaluated using standard expression ]lOt . 

The emission of neutrons during the temporal evolution of the trajectory is simulated as follows. At each time step, 
the probability of emission of a neutron, T/r„ (where t„(= ?i/r„), r are the neutron decay time and the time step 
of the calculation, respectively), is computed and compared with a random number i?jv from a uniform probability 
distribution. The emission of a neutron is assumed to take place, if it satisfies the following criterion pO| ; 

T/r„ > Rn. (3.10) 



If the condition (3.1C) is not satisfied, no emission of neutron takes place. The time step r is chosen in such a way 
that it satisfies the condition r/r„ <^ 1. Consequently the probability of emission of two or more neutrons in time 
r would be extremely small. The calculation is continued over the whole trajectory for a number of times at each 
angular momentum £ to estimate the average prescission multiplicity at each £. The calculation is then repeated for 
all allowed values of angular momentum to compute the aveage value of prescission neutron multiplicity Upre [ p^ . 

The calculated values of Up^e have been displayed in Fig. ^ alongwith the respective experimental data as a function 
of the initial excitation energy of the compound nucleus for two different mass regions. The solid curves are the results 
of the present calculations and the symbols correspond to experimental data ||l8|-po|| . It is seen that for heavier systems 
{AcN ~ 200) (lower half), the theoretical predictions are in good agreement with the corresponding experimental 
data. For lighter systems {Acn ^ 150) (upper half), the experimental values of ripre have larger uncertainties and 
fluctuations, and the theory is seen to reproduce quite well the average trend of the data. A part of this fluctuation 
in neutron emission here may be due to specific structure effects of different compound systems; for example, ^^^Yb 
(filled diamond) is quite neutron deficient compared to ^^^Yb (open triangle), and neutron emission from the former 
is therefore expected to be somewhat less. Similarly, at high incident energies (> 10 MeV/nucleon) the observed 
multiplicity (open diamond) was found to be lower than the average theoretical trend, which may be due to the 
noninclusion of the effect of preequilibrium emission in the present calculation. 

The fragment mass asymmetry dependence of neutron multiplicity is displayed in Fig.^jfor ^^O (Eiab = 158.8 MeV) 
induced reactions on ^^"^Sm, ^^'^ Au and "^^^JJ ||l^. The solid circles correspond to the experimental data and the solid 
lines are the theoretical predictions of the same. It is found that the present calculations agree quite well with the 



experimental data in all the cases. The value of the constant K (Eqn. 3.2) was found to be 161±3 which is independent 
of the mass of the compound sys tem. It is, therefore, interesting to note that with the inclusion of the term h(a) in 
the friction form factor (Eqn. |3.2[) , we are able to explain the prescission neutron multiplicity data for both symmetric 
as well as asymmetric fission with the same value of the viscosity coefficient, (= 4 x 10~^^Mey • sec ■ fm^^). 



2. Energy of emitted neutrons 

The kinetic energy of the emitted neutron is extracted through random sampling technique pO| . Assuming that the 
system is in thermal equilibrium at each instant of time t, the energy distribution of the emitted neutrons is represented 
by a normalised Boltzmann distribution corresponding to the instantaneous temperature of the system. From a 
uniformly distributed random number sequence {a;„} in the interval [0,1], another random number sequence {j/n} 
with probability distribution f(y) is constructed, where f(y) ^ exp(—f3(t)y) is a normalised Boltzmann distribution 
corresponding to the temperature (3(t) at any instant of time t. Then, the sequence {?/«} is obtained from the sequence 
{x„} by the relation, 

y(x)=F-Hx). (3.11) 

Here, is the inverse of the function F(y) = x — f(y)dy, which is computed numerically by forming a table of 
integral values. The energy of the emitted neutron is given by En — yE^"-^ , where i?™"^ is chosen in such a way that 



7 



the Boltzmann probability at that energy is negHgible for all instants of time t. After the emission of the neutron, 
the intrinsic excitation energy is recalculated and the trajectory is continued. 

The average energy of the prescission neutrons, < En > has been plotted as a function of the compound nuclear 
mass, AcN, in Fig. It is seen from the figure that the theoretical predictions of < £"„ > [solid curve) are in good 
agreement with the respective experimental data (filled circles). 



C. Average and variance of TKE 



The temporal evolutions of the variables g{t), C{t), G{t) along the fission trajectory have been computed for a 
representative system ^^O + ^^''Sn and the results are plotted in Fig. |[ It is seen from the figure that, initially, all 
of them increase steeply and then their magnitudes become nearly constant throughout the rest of the trajectory. 
Furthermore, the calculation shows that C^{t)/g{t)G{t) <C 1, which implies that the correlation of position and 
velocity of the elongation variable (r) is much smaller compared to their respective variances. The variance of energy 
and average of total kinetic energy (TKE) at scission point are given by, 

4 = (M0'.9(i) + [d{Vc + VN)/dr]^G{t), (3.12a) 
< E{t) > = fig{Uc)/2 + Edet- (3.12b) 

The contribution of term (typical ly ~ 2 iir(d(Vc + VN)ldr)C{t) ) involving the correlation between position and 
velocity has been neglected in Eq. ( 3.12a[ ) as it is quite small compared to the other terms invoving the variances 



of position and velocity. The quantity tsc is the time at scission point and Edet is the deterministic value of total 
fragment kinetic energy (TKE) after scission and ~ 100 - 200 MeV. It is assumed that the variation of the potential 
over the narrow width of the probability distribution is small so that the average of the potential is approximated as 
the value of the potential at the mean position. The variation of the kinetic energy variance (t|; as a function of time 
has also been displayed in Fig. IJ. The value of a^^ is also seen to increase steeply at the beginning and then it becomes 
nearly constant throughout the rest of the time. As envisaged earlier, the result clearly shows that cje/ < E > <^ 1, 
which demonstrates the validity of asymptotic expansion in deriving the result instead of solving the Fokker-Planck 
equation in detail. 

The theoretical predictions of aE{th) and mean total kinetic energy (TKE) (solid curves) for the fission of several 
compound systems produced in the 158.8 MeV ^^O, 288 MeV ^^O induced reactions on various targets have been 
displayed in Figs. |a, respectively, alongwith the experimental data |l^ (filled circles). From the figure it is 
observed that theoretical predictions of TKE agree quite well with the respective experimental data for all the 
systems studied. However, it may be noted here, that the experimental values of asiexp) are usually obtained by 
averaging over the full mass yield spectrum. Therefore, asiexp) consists of two terms, viz., (i) contributions arising 
due to stochastic fluctuations in the dynamics of fission process, ctej and (ii) contributions from the variation of the 



mean kinetic energy with the fragment mass asymmetry, (jE{kin). So, asiexp) may be written as |21|, 

a%{exp) = a'^ + a'^{kin), (3.13a) 
oI^Y.''e{Ai,A2)-Y{A,), (3.13b) 

alikin) = ^[^- < £^(^1,^2) >]2 • Y{Ai). (3.13c) 

Ai 

Here a%{Ai, A2) and < E{Ai,A2) > are the variances and mean values of the total kinetic energy of two fission 
fragments with mass numbers Ai and A2 (compound nucleus mass Acn = Ai + A2), E being the average of < 
E{Ai,A2) > over the normalised fragment mass yield, Y{Ai) with J^Ai^i^-i) = 1. Thus, the calculated value of 
a'^(th) is then compared with the stochastic component of the experimental variance, i.e., which is obtained after 
substracting a%(kin) from (j\{exp) (as mentioned above). 

We have extracted a'^{kin) for a few systems for which the experimental fragment mass yield data are available 
[Q, taking < E{Ai,A2) > from Viola systematics [|2|. The values as, i-e., a/ {a'^{exp) — a^{kin)) , are shown in 
Fig.^ as open triangles and they agree very well with the predicted values of TKE variance. It is seen from Fig. |^ 
that when the projectile energy (and vis-a-vis the excitation energy of the fused composite) is relatively lower, the 
calculated values are in fair agreement with the data. However, the calculation underpredicts the experimental value 
of aE for the heaviest target considered (^'^^U in the present case). With the increase in the projectile energy (and 
the excitation energy of the composite), the theoretical predictions are found to underestimate the corresponding 
experimental values and the discrepancy between the two increases with the increase in mass number (Fig. pb),. 
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We have also studied the fragment mass asymmetry dependence of energy variance, <t'^{Ai, A2) for some represen- 
tative systems and the results are displayed in Fig. g. It is seen from Fig. g that the theoretical values of variances 
have only a weak dependence on the fragment mass asymmetry. 



IV. SUMMARY AND CONCLUSIONS 



We have developed a generalised formulation of asymptotic expansion of the Fokker-Planck equation for the systems 
where the diffusion coefficient depends on the stochastic variable explicitly. With the assumption that the relative 
fluctuation of collective variable is small wc have derived the equation for various moments. The formalism is applied 
to the case of fission where the fluctuation in total kinetic energy is small as compared to its mean value. We have 
taken only one degree of freedom, namely the elongation axis in our calculation. However, one could incorporate neck 
degree of freedom also in a more realistic calculation based on the present formalism. The primary motivation of the 
present work is to show that this formalism could explain the basic features of the fission dynamics quite satisfactorily 
without invoking the solution of Fokker- Planck equation or Langevin equation in detail. 

The present model is found to explain fairly well the observed neutron multiplicities and their fragment mass 
asymmetry dependence as well as the average energy of the evaporated neutrons over a wide range of mass and 
excitation energies of the compound system with a single value of the viscosity coefficient, fiQ. The predicted values 
of TKE are found to be in good agreement with the experimental data and the theoretical estimates of the associated 
TKE variances are also found to agree quite well with the respective numbers extracted from the experimental data 
for the systems where the fragment mass yield data are available. For a more direct test of theoretical models it is 
necessary that experimental estimation of variances should not have admixture of other contributions arising due to 
the variation of mean kinetic energy over different mass yields. This may be achieved if measurements are done in 
smaller mass bins. 

In the present studies, the correlation of the position and velocity of the elongation axis has been found to be small. 
However, in the cases where such condition is not valid the energy variance still can be calculated by adding a term 
2fi'r{d{Vc + Vj\r)/dr)C{t). The procedure developed here could systematically generate higher order hierarchies for 
relatively larger fluctuations than the ones encountered in the present studies. In those cases one may have to solve 
the higher order equation which would involve higher order derivatives of the functions h(x, y) a,ndH{x, y), in general. 
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FIG. 1. Prescission neutron multiplicities plotted as a function of the initial excitation energy E*„i of the compound nuclei 
of masses Acn ~ 150 {upper half), and Acjv ~ 200 (lower half). The solid curve is the present calculation. Different symbols 
correspond to different sets of experimental data, (ie^fiUed circle |^], open inverted triangle filled inverted triangle [ |l8[ , 
open triangle [Q, filled triangle pof , open diamond [^, filled diamond jTsj, open square [^, filled square ]l^). 
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2. Prescission neutron multiplicity Upre as a function of fragment mass asymmetry aasyior ^^O induced reactions on 
^^'^ Au and ^^*U. Filled circles correpond to the experimental data |l^ and the solid curves are the present calculations. 
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FIG. 3. mean energy of the evaporated neutrons (< En >) plotted as a function of Acn- The sohd curves are the present 
calculations, and the filled circles are the corresponding data |18|. 
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FIG. 4. Variation of (a) g{t), (b) C{t), (c) G{t) and (d) a% as a function of time t for the system ^'^O+^^^Sn. 
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FIG. 5. Variation of as and mean TKE as function of target mass number At, for (a) 158.8 MeV , and (h) 288 MeV 
^^O induced fission reactions. Filled circles correspond to the experimental data |l8| and solid curves are the present theoretical 
results. Open triangles are the values of ctb obtained using Eqs. (3.13) (see text). 
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FIG. 6. Variation of predicted values of ctb as a function of fragment mass asymmetry, asym ~ \Ai — A2\/Ac 
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